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ABSTRACT 

We present results from a study of molecular gas and dust in the vicinity of the 
Herbig Haro object HH2. Emission from the sub-mm continuum, 12 CO and HCO + was 
mapped with angular resolutions ranging from 14 " to 5 " (or 0.01 pc at the distance of 
HH2). The continuum shows an extended dust clump of mass 3.8M Q and temperature 
22K, located downstream of the bright optical HH knots. However, a compact emission 
peak lies within 0.01 pc of the low-excitation H 2 -promincnt shocks, with a luminosity 
consistent with local heating by the outflow. 

The HCO + emission shows two velocity components: firstly, ambient-velocity gas 
lying in a region roughly corresponding to the dust clump, with abundance enhanced 
by a factor of a few close to the H2-prominent knots. Secondly a component of high- 
velocity emission (20 kms -1 linewidth), found mainly in a collimated jet linking the 
low-excitation HH objects. In this high- velocity jet, the line wings show an abundance 
ratio Xhco+ Ixco oc v 2 , with an HCO + enhancement compared with ambient gas of 
up to ~ 10 3 at the most extreme velocities. Such high abundances are consistent with 
models of shock chemistry in turbulent mixing layers at the interaction boundaries of 
jets. Extrapolating this effect to low velocities suggests that the more modest HCO + 
enhancement in the clump gas could be caused by low velocity shocks. A UV precursor 
may not, therefore be necessary to explain the elevated HCO + abundance in this gas. 
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1 INTRODUCTION 

The Herbig-Haro objects HH1 and HH2 are two of the 
brightest optical shocks associated with outflows from young 
stars. Separated by 140 " or 0.3 pc at the distance of 460 pc, 
both their proper motions and the detection of a faint jet in- 
dicate they are driven by the embedded 32 Lq Class I young 
star known as VLA1 (Herbig & Jones, 1981; Pravdo et al., 
1985). Because of their brightness, HH1 and 2 have been the 
subject of intense study, particularly in optical and near- 
infrared lines used to tracing shocked gas (eg Hester et al., 
1998). The images show complex and sometimes bewildering 
structures which, in the case of HH1, are often interpreted as 
multiple overlapping bow shocks. HH2, however, shows an 
apparently more random set of knots, although it too has a 
number of mini-bow shocks identified through infrared im- 
ages (eg Davis et al., 1994). In order to aid orientation in 
the region, in Figure 1 we show a near-IR H2 image of the 
cluster of shock fronts that comprise HH2 (from Davis et al. 
1994). Individual features and knots are labelled. 

In addition to HH1 and 2, a pair of much larger, 



more diffuse HH objects HH401/402 can be seen almost 
20' (2.6 pc) further from VLA1 at the same position an- 
gle (Ogura, 1995). Thus HH1 and 2 may actually be inner 
knots in a "parcsec-scale" outflow. 

Single-dish HCO + observations in the immediate vicin- 
ity of HH2 have shown an emission peak downwind of the 
optical knots, and it was suggested that this could be due 
to enhancement in the HCO + abundance in the ambient gas 
(Davis et al., 1990; Dent 1997). Similar HCO+ peaks have 
been found in several other outflows (eg Rudolph & Welch, 
1992). A common feature is that the emission lies close to 
ambient velocity, and the peak appears ahead of the HH 
shocks. It was proposed that these clumps are caused by 
UV photons from the HH shock inducing non-equilibrium 
chemistry in ambient material in either purely the gas phase 
(eg Wolfire & Konigl, 1993), or after the release of icey 
grain mantles (eg Taylor & Williams, 1996; Viti & Williams, 
1999). It is also possible for HCO + in entrained shocked gas 
to be chemically enhanced by the shock itself (eg Taylor & 
Raga, 1995). Further imaging of HH1 and 2 showed compact 
clumps of NH3 in the vicinity (Torrelles et al., 1992), again 
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Figure 1. Image of H2 emission from the vicinity of HH2, based 
on data from Davis et al. (1994). The individually-identified knots 
are labelled and the flow direction is indicated by the dashed 
arrow in the upper right. Note that "Knot 8" is only seen in H2 
emission (it has no known optical counterpart). 

thought to be due to changes in the molecular abundance. 
A more extensive survey of other species towards the HCO + 
peak has shown abundances consistent with the effect of UV 
photons from the nearby HH shock (Girart et al., 2002). 

The molecular outflow associated with HH1-2 is not 
prominent, but line wings have been detected in CO to- 
wards HH2, albeit at relatively low velocities (eg Correia et 
al., 1997; Dent, 1997). Maps of the line wing emission over 
the whole region show a bipolar outflow at certain velocity 
ranges, which appears to terminate near HH1 and HH2; it 
was suggested that the flow axis is inclined only 5-10°from 
the plane of the sky (Moro-Martin et al., 1999). 

In order to investigate the relationship between the am- 
bient clump and high-velocity gas, we have mapped the re- 
gion in the sub-mm continuum and J=3-2 CO line, and also 
obtained higher-resolution J=l-0 HCO + images. 



2 OBSERVATIONS 

The continuum observations were made using SCUBA on 
the JCMT in June 1998. Imaging was carried out at 850 fim 
using a standard 64-position fully-sampled jiggle map, and 
data were calibrated using Mars, resulting in a calibration 
accuracy of ~ 10%. A fully-sampled map of the J=3-2 12 CO 
line was also obtained using the JCMT; additional J=3-2 
spectra of the isotopomers 13 CO and C 18 were extracted 
from the JCMT Archive. These data were taken in 1994, 
using the receiver RxB3i and the DAS spectrometer with a 
spectral resolution of 0.3fcms -1 . The beamsize of the JCMT 
in both these continuum and spectral line observations was 
14 ", and estimated pointing uncertainties were less than 2 ". 

We conducted interferometric imaging of HCO + J=l-0 
line emission (rest frequency 89.18852 GHz) using the 6- 
element Nobeyama Millimeter Array (NMA) in 1999 March 



and May. The primary beam size (field of view) is 83 " 
(fwhm) at 89 GHz. Observations were performed with the 
C and D configurations, and the resulting synthesized beam 
size was 6.9" x 4.6" at P.A.= —20°. Since the minimum 
projected baseline length was 3.9 kA, our observations were 
insensitive to structure extended by more than ~ 52 . The 
phase tracking center was set on the position of the peak of 
J=3-2 HCO+ (5 h 36 m 26.91 s , -6° 47' 31.9"; 2000.0). We 
employed SIS receivers which had double sideband system 
noise temperatures of 200-300 K toward zenith. We used 
the high spectral resolution FX correlator for the HCO + 
emission and the Ultra- Wideband Correlator for calibration 
using continuum emission. The FX correlator provides a ve- 
locity resolution of 0.1 km s _1 and total velocity coverage of 
108 km s -1 at this frequency. We used 3C273 as the band- 
pass calibrator and 0528+134 as a phase and gain calibrator. 
By comparing with Uranus through 3C84 and 3C454.3, we 
estimated the flux density of 0528+134 to be 2.5Jy during 
both observing periods, with an estimated 10% uncertainty. 
Imaging was performed using the AIPS package to make 
clean maps. The typical resulting rms noise levels were 75 
mjy beam -1 with a resolution of O.lfcms -1 . 



3 RESULTS 

3.1 Dust continuum 

The emission from HH2 observed at 850 /im is shown by 
the contour map in Figure 2. For comparison purposes, the 
shocked H2 image from Davis et al. (1994) is superimposed 
as a greyscale. These H2 knots lie towards the end of the 
Southeast lobe of the outflow; the flow origin (VLA 1) is 
situated some ~ 90" to the Northwest. An extended region 
of sub-mm continuum emission, of dimensions 60 " x 40 " 
(0.13x0.09 pc, fwhm) is clearly seen downstream from most 
of the HH objects. The peak in sub-mm emission lies at 
5 h 36 m 27.0 s ,-6°47'21" (2000.0), and could be either due 
to a local maximum in the dust temperature or density. The 
J=4-3 HCO + line observed with the same spatial resolu- 
tion peaks 10" to the south (Dent, 1997), and so the posi- 
tions are consistent within the estimated errors. Girart et al. 
(2002) also found that the SO 32 — 2i line reaches a peak at 
this position. Although there are at least 3 bright sub-mm 
continuum clumps within 2 arcmin of VLA 1 (Chini et al., 
2001), most of these have evidence of an embedded young 
star. The dust clump near HH2 has no such known object, 
even though the extinction through the cloud is relatively 
low (see Girart et al., 2002). The close association with the 
optical, infrared and HCO + emission suggests that the sub- 
mm continuum is tracing a cloud closely related to the HH 
objects. Furthermore Fig. 2 shows that at least two faint HH 
knots are apparently embedded within the sub-mm cloud; 
the significance of this will be discussed later. Also it is clear 
that the cloud has a much sharper edge facing VLA1, sug- 
gesting that the stellar jet has eroded away one side of the 
cloud. 

The total integrated 850 /im flux from the extended 
cloud associated with HH2 is estimated at 2.1+0.4 Jy (this 
is measured over a region of 80 " x 60 arcsec at PA = 30° 
centred on the peak). The peak flux is 0.20+0.04 Jy, con- 
sistent with an earlier single point photometric observation 
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Figure 2. Contours of sub-mm continuum emission in the vicin- 
ity of HH2, superimposed on the near-infrared H2 image from 
Davis et al. (1994). Contour interval is 22 mjy/beam, with low- 
est contour 44 and peak of 200 mjy/beam. The beam size in the 
sub-mm data is 14" fwhm. VLA1, the driving source for the flow, 
is shown by the cross in the upper right; the outer edge of the 
dust disc associated with this young star can be seen towards the 
upper right of the contour map. The map coordinates arc J2000.0. 
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Figure 3. Spectral energy distribution of the cloud near HH2. 
The 850 (im point represents the total emission from the cloud 
associated with HH2 (see text). Data points indicated by stars 
are based on results from Molinari & Noriego-Crespo (2002), and 
open circles are IRAS data points taken from Pravdo & Chester 
(1987). The solid line shows a fit to the long-wavelength data 
only, with T=22K, /3=1.5. The 12 and 25[im points are fitted 
assuming the same value of /3, and with T=220K; however, this 
mid-infrared emission may arise from a nearby source not directly 
associated with HH2 - sec text for details. 



at 800 /im (Dent, 1997). By comparison, the exciting star 
VLA 1 has a peak flux of 3.3 Jy (Chini et al., 2001). Both 
HH2 and VLA 1 show a compact central source superim- 
posed on an extended low-level plateau; in both cases the 
latter may be heated by the ambient interstellar radiation 
field rather than the central source (see Chini et al., 2001). 

Pravdo & Chester (1987) detected emission in the re- 
gion of HH2 at 12 and 25 /im using IRAS, which they ten- 
tatively ascribed to the HH object. However, their maps 
indicated the mid-infrared peak may be ~ 1 ' west of HH2; 
this is supported by more recent ISO data (Cernicharo et 
al., 1999), suggesting the mid-infrared continuum is conci- 
dent with the H a rim noted by Reipurth et al. (1993) to lie 
50 " west of HH2. However, far-infrared ISO results did show 
continuum emission from the region of HH2 itself (Molinari 
& Noriega-Crespo, 2002); the estimated contamination from 
the bright source VLA 1 was <10%. Figure 3 shows the 
integrated continuum fluxes. A grey body fit to the long- 
wavelength data gives a temperature of 22±2 K and dust 
opacity index, /3=1.5±0.2; this would imply a total cloud lu- 
minosity of 13±5 Lq. The mid-IR points are consistent with 
T=220K, but these are ignored in this and subsequent anal- 
ysis for reasons given above. Furthermore the temperature 
of this component is sufficiently high that emission is un- 
likely to significantly affect the far-infrared fluxes and hence 
the fit (see Fig. 3). However, further high-resolution imaging 
in both the far and mid-infrared is required to confirm this 
result. 

The NH3 observations of Torrelles et al. (1992) indi- 



cated kinetic temperatures in the region of the HH2 clump 
of <20K. Girart et al. (2002) derive CO excitation temper- 
atures of ~13K. The peak brightness temperature in 13 CO 
and 12 CO is ~10K (see below), which is a lower limit to the 
kinetic temperature assuming the line to be optically thick 
and in LTE. 

If we adopt a temperature of 22K, the total mass of 
the HH2 cloud is 3.8±0.4Mq, assuming a gas:dust ratio 
of 100, and an 850 /im dust mass opacity of W~ i m 2 kg~ 1 
(eg Henning et al., 1995). This gives a mean gas column 
density of 1.8 ± 0.2 x 10 22 cm -2 and mean space density of 
5 ± 0.5 x 10 4 c?n~ 3 (assuming a symmetrical cloud). Note 
that the errors on these derived parameters may be larger - 
perhaps up to a factor of 2 - because of the uncertainty in 
the mm- wavelength dust emissivity (eg Henning et al.). 

The contours in Figure 2 show a clear emission peak 
near HH2; assuming the temperature and dust opacity here 
is the same as that of the whole cloud, then the luminosity 
of this component is estimated as ~ 1.2 ± 0.3 Lq. Its mass 
is 0.3 ± 0.08 M , peak column density 2.3 ± 0.6 x 10 22 cm~ 2 
and peak space density is 2 ± 0.5 x 10 5 cm -3 . These den- 
sities are very similar to those found from recent observa- 
tions of several molecular species towards HH2 (Girart et 
al., 2002), although the derived temperature is somewhat 
higher. This general agreement confirms that the dust and 
molecular gas lie in the same warm clump, and also sug- 
gests that the sub-mm dust emission characteristics are not 
significantly affected by the nearby shock. 



© 0000 RAS, MNRAS 000, 000-000 



4 W. R. F. Dent, R. S. Furuya, C. J. Davis 




Figure 4. Spectra of CO J=3-2 towards the continuum peak 
near HH2. From top to bottom, these are 12 CO, 13 CO and C ls O. 
Intensity scale is T m (,, and velocity (in fcms -1 ) is with respect 
to Isr. The vertical line represents the systemic radial velocity of 
the molecular clump associated with HH2. 

3.2 Single-dish CO 

The outflow axis from VLA1 lies close to the plane of the 
sky, resulting in relatively faint CO line wings. Nevertheless 
a collimated molecular outflow has been observed close to 
VLA1 (Correia et al., 1997; Choi & Zhou, 1997), and fainter 
red-shifted CO extends as far as HH2 (Moro-Martin et al., 
1999). In an attempt to understand this termination region, 
we observed HH2 itself using the J=3-2 transition of CO 
and its isotopomers; the higher energy level of this transi- 
tion compared with earlier observations means it is more 
sensitive to the warmer gas. Figure 4 compares spectra at 
the location of the continuum dust peak. The C ls O line can 
be fit by a single Gaussian component at 7.25 ± 0.2fcms _1 , 
close to that of the HCO+ (7.0fcms _1 ; Davis et al., 1990), 
CI (Dent, 1997), and other species such as SO (Girart et al., 
2002). This velocity is significantly different from the gas 
around VLA1, which has a radial velocity of 10.5 kms^ 1 
(Choi & Zhou 1997). 

The ratio of peak brightness temperature of 13 CO to 
C 18 is ~5 which, assuming a similar emitting area and LTE 
conditions, indicates that the C 18 line is optically thin. The 
integrated C 18 intensity is 2.6Kkms~ 1 ; with an excitation 
temperature of 22K (see above) and a C ls O abundance of 
10~ 7 , this gives a total gas column density of 1 .4 x 10 cm 
This is consistent with the value estimated above from the 
sub-mm dust emission. 

The 12 CO line in Figure 4 shows distinct high- velocity 
wings at relative velocities |« re ;| > 3fcms _1 . The spatial dis- 
tribution of this high-velocity gas in the vicinity of HH2 is 
illustrated in Figure 5; both the red and blue-shifted emis- 




tours delineate red-shifted gas integrated over the range 13 - 18 
kms^ 1 and the solid contours show gas in the range 0-5 kms^ 1 . 
Contour interval is 0.8 Kkms~ 1 , with a base of 0.8 (temperature 
scale in antenna units). The coordinates are J2000.0; beam size 
is 14" (fwhm). The location of some of the major HH clumps arc 
indicated by triangles; from the South these are HH2L, Knot 8, 
HH2F, E and A. The peak of the sub-mm continuum is shown by 
a cross and the presumed direction of the outflow from VLA1 is 
shown by the arrow. 

sion peak within 10 " of the brightest Hb knot HH2E. A 
second component of blue-shifted gas lies south of HH2L. 
It was shown by Moro-Martin et al. that red-shifted CO 
around HH1-2 traces a large collimated flow centred on 
VLA1; high-resolution observations of blue-shifted HCO + 
(to be described below) show that it too lies in a collimated 
jet pointing back towards VLA1. It is thought unlikely that 
the dust clump in Fig. 2 harbours a young outflow source, 
as the low extinction through the cloud would render it de- 
tectable in the infrared. So the abrupt termination of the 
outflow near HH2E (Fig. 5), suggests this is the dominant 
interaction region of the jet and ambient clump. However 
the spatial resolution of these single-dish data is insufficient 
to compare accurately with the infrared image. 

3.3 Interferometric HCO+ 

The HCO + J=l-0 emission from the vicinity of HH2 is 
bright and complex, but the spectra have two distinct ve- 
locity components. Firstly, bright and relatively narrow 
ambient-velocity emission with a profile similar to that of 
13 CO; this is illustrated in Figure 6a which shows the line 
towards the continuum peak. Secondly a high-velocity wing 
component, with a full width of 20 kms^ 1 (to the noise 
level). This dominates in the region of HH2E (see Figure 6b). 

Both HCO + spectra show an absorption dip at vi sr ~ 
10/cms -1 , similar to the velocity of the molecular gas near 
VLA 1 and in much of the Orion region (Davis et al., 1990; 
Choi & Zhou, 1997). It is likely that a cool region of this 
cloud along the line of sight is absorbing the emission from 
HH2. However, a second dip in the HCO + spectrum at 
vi S r ~ 7fcms~ 1 is likely to be cool foreground gas associ- 
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Figure 6. Comparison of the J=l-0 HCO~*~ line (thick his- 
tograms) with J=3-2 12 CO (dashed lines) towards (a) the contin- 
uum peak and (b) HH2E. High-velocity emission is seen in both 
molecules, although the HCO+ line is wider than CO. The HCO+ 
data have been smoothed to the same resolution as the single-dish 
data (14"); intensity scale is T m (, and velocity frame is Isr. 



ated with the HH2 clump itself, as seen in the CO spectra 
in Fig. 4. 

The map of ambient velocity HCO + (Figure 7) shows 
a main peak of size 20" (fwhm), consistent with the lower- 
resolution observations of Davis et al. (1990). Fig. 7 com- 
pares the data with the sub-mm continuum emission and the 
locations of some of the prominent HH knots (c.f. Fig. 1). 
The gas and dust distributions are similar but not identi- 
cal; notably the HCO + shows two compact peaks within 5 " 
(O.Olpc) of HH2L and Knot 8, whereas the dust is seen in a 
broader North-South ridge extending over ~ 50". It is pos- 
sible that some of the large-scale HCO + emission (> 52") 
has been resolved out by the lack of short-baseline data from 
the interferometer. Additional instrumental effects due to in- 
complete UV coverage can be seen as negative contours in 
Figure 7, although these are at a relatively low level com- 
pared with the peak intensity. The map may also be affected 
by the self-absorption near the line centre, although this is 
relatively narrow compared with the main emission profile. 
Despite these potential difficulties, it does appear that the 
ambient- velocity HCO + emission closely traces the dust con- 
tinuum; however, it appears relatively brighter within O.Olpc 
of the HH shocks HH2L and Knot 8. This enhancement will 
be discussed further in Section 4.1. 
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Figure 7. Contours of ambient-velocity HCO+ (5.0 to 8.5 
kms^ 1 ) superimposed on a greyscale image of the sub-mm con- 
tinuum. Contour levels are -8, -4, 4, 8, 12 ... Kkms -1 ; negative 
contours are shown in black. Note that the negative contours 
are due to the limited UV coverage of the observations; they are 
not thought to significantly affect the morphology of the peak 
emission. Also shown by triangles are the locations of the HH 
objects HH2L, Knot 8, HH2F, HH2E, and HH2A (in order, from 
the South). Coordinates are J2000. The synthesised beam size is 
6.9" x 4.6" at PA = -20°. 



3.4 High-velocity HCO+ 

Maps of the integrated high-velocity HCO + emission arc 
shown in Figure 8, with the near-infrared H2 image (Davis 
et al., 1994) superimposed. The distribution of high- velocity 
gas appears similar to that of CO shown above, after ac- 
counting for the differing resolution and relative insensitivity 
of the interferometric observations to the large scale struc- 
ture (> 52"). High-velocity red and blue-shifted HCO+ is 
found within 5 " of the HH objects HH2K, E, F and also 
Knot 8. Notably, these knots are relatively bright in shocked 
H2 compared with optical lines such as H a , whereas HH2A 
and HH2H, which dominate optical images, are relatively 
faint in H2 (eg Noriego-Crespo & Garnavich, 1994). Multi- 
level line analysis shows that they also have lower excitation 
temperatures (Eisloffel et al., 2000). There is weak evidence 
of high- velocity red-shifted gas near the western HH clump 
HH2A. 

The images of high-velocity HCO + (particularly the 
blue-shifted gas) show that the northern section of HH2, 
rather than being a jumble of shocks, actually forms a rather 
coherent well-collimated jet. This includes HH2K, E, and F, 
terminating in Knot 8 (and thus confirming that the latter 
object, which is only detected in H2, is part of the flow). The 
Position Angle of this jet measured from VLA1 is 139° - sig- 
nificantly different from the PA of the optical peak HH2A 
measured from VLA1 (150°). By comparison, the apex of 
the large-scale outer bow shock HH402 lies at PA w 138° 
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(Ogura ct al., 1995). Assuming a constant flow direction, 
these results suggest that the dominant outflow position an- 
gle is 138°. This is further confirmed by the fact that the 
sub- mm continuum peaks at PA ~ 140° from VLA1; if we 
assume it is heated by the jet, this traces where most of the 
outflow energy is being deposited. 

Giarart et al. (1999) suggested that HCO + may be used 
to identify shocked gas in sections of outflows that are no 
longer visible in optical or infrared lines. Note that the cool- 
ing times associated with the optical/IR lines are very short 
- of the order of a few years - so these emission lines can 
rapidly fade. The HCO + lines, on the other hand, emit from 
molecular material that is already relatively cool, where 
steady depletion via dissociative recombination with elec- 
trons may take > 10 2 years (Nejad & Wagenblast 1999), 
although this is dependent on the electron abundance. In- 
deed, the apparent HCO + abundance enhancement evident 
in our data may be limited by the lifetime of the clump in 
HH2 which, if we assume is destroyed by the passage of the 
shock front, will be of the order of a few thousand years. 



4 DISCUSSION 

4.1 The ambient clump 

Although the sub-mm results show an extended cloud down- 
wind of the bright optical HH objects, the central sub-mm 
peak coincides with a number of low-excitation HH knots. 
The cloud has no evidence of an internal young star, so how 
is the dust heated? Either it could be through the inter- 
stellar radiation field in the young cluster, by UV radiation 
from optically-bright knots such as HH2A, or by shock heat- 
ing. The fact that the submm continuum peaks so close to 
HH2 suggests the heating - at least of this central peak - 
is local, rather than from the external interstellar radiation 
field. One possible local energy source is UV from the HH 
shocks, as used to explain the HCO + overabundance in the 
clump (eg Wolfire & Konigl, 1993; Taylor & Williams, 1996; 
Viti & Williams, 1999). Bohm-Vitense et al. (1982) found 
the dominant UV source in HH2 lies within ~ 6 " of the op- 
tically bright high-excitation knots HH2H and A. Thus the 
main UV source is ~ 25 " Northwest of the sub-mm peak, 
and the fraction of UV radiation intercepted by the dust 
clump would be ~ 5%. The total de-reddened UV luminos- 
ity of HH2A, estimated from Fig. 2 of Bohm-Vitense et al. is 
~ 0.3 L©, which would result in an intercepted luminosity at 
the sub-mm peak of ~ 0.015 L©, considerably less than that 
observed (~ 1.2 L©). Clearly the UV flux from HH2H/A is 
inadequate to heat the dust peak. 

Alternatively the clump may be warmed by local shocks 
from the outflow jet. Moro-Martin et al. give the total out- 
flowing mass from VLA1 as 0.2 M©; with a mean flow ve- 
locity of 30fcms _1 (deprojected assuming an inclination of 
10°) and age of W 4 yrs, this results in a total flow mechan- 
ical luminosity of ~ 1.0 L©. This is adequate to explain the 
local heating of the sub-mm peak, although not of the whole 
cloud shown in Fig. 2. 

Both the sub-mm continuum and the ambient-velocity 
HCO + reach a maximum within 5 " (0.01 pc) of HH2L and 
Knot 8. To estimate the HCO + enhancement in this region, 
we can compare the molecular abundance with a potentially 
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Figure 8. Contours of high-velocity red (upper panel) and blue- 
shifted HCO+ superimposed on the H2 image. Red-shifted gas 
is integrated over +11.0 to +20.0 krns^ 1 and blue-shifted gas 
from -5 to +5.0 ferns" 1 . Contours arc -0.4, -0.2, 0.2, 0.4, 0.6... 
Kkms^ 1 . Map coordinates arc J2000. The synthesised beam size 
is 6.9" x 4.6" at PA = -20°. 



more "benign" region of the cloud 25 " to the South. Data 
are smoothed to the same spatial resolution (14" fwhm), 
and we make an initial assumption of an excitation tem- 
perature of 22K (see above), LTE conditions and optically 
thin emission. Towards the dust peak and southern region, 
the HCO+ column densities (N L (HCO + )) are 3.6 x 10 13 
and 1.1 x 10 13 cm~ 2 respectively. This compares with esti- 
mates derived from the dust of Nl{H2) = 2.3 x 10 22 and 
1.8 x 10 22 , giving abundances, Xhco+ °f 1-6 x 10~ 9 and 
6.5 x 10~ 10 . Girart et al. (2002) used multi-transition mod- 
elling to show that the HCO + excitation temperature may 
be somewhat less than 22K, possibly because the mean cloud 
density is less than the critical density of 10 6 cm -3 . Further- 
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more they suggest the line may be optically thick towards 
the bright peak, which is supported by the self-absorbed 
profile in Fig. 6a. But if we assume similar conditions ap- 
ply in the two regions, it indicates Xhco+ m the ambient 
gas is a factor of >2.5 higher within ~ 0.01 pc of HH2L, 
compared with that 0.05 pc to the South. Furthermore, the 
abundance in the Southern region is similar to that in many 
other quiescent molecular clouds (eg Nejad & Wagenblast, 
1999). 



4.2 High-velocity molecular gas 

High-velocity HCO + and CO emission is closely associated 
with a line of H2 knots terminating at Knot 8 (Fig. 8). This 
suggests a single coherent jet impacting the dust clump, 
rather than a broad wind with multiple shocks as might 
be suggested by the complex optical images. A similar spa- 
tial association of high-velocity HCO + with H2 peaks was 
found in NGC2071 (Girart et al., 1999) and DR21 (Garden 
& Carlstrom, 1992), and it was proposed that the HCO + 
enhancement is due to ion-chemistry reactions in the low- 
velocity C-shocks. 

After accounting for the differing resolution of the data, 
there appears to be no significant difference in the spatial 
distribution of high- velocity CO and HCO + . However, there 
is a clear difference in the line brightness ratio as a function 
of velocity, illustrated in Figure 9. We find this ratio changes 
from ~0.1 at the lowest relative velocities likely to be un- 
contaminated by ambient emission (\v re l\ ~ 3fcms _1 ), up 
to ~2-3 at the highest velocities (\v re l\ ~ 12fcms _1 ). Here 
we use v re i = vi sr — v sys , where v sys , the systemic clump 
velocity, is taken to be 7.Qkms~ . Note that the inclina- 
tion of the HH1-2 flow is only ~ 10° to the plane of the 
sky (eg Moro-Martin et al., 1999), so the absolute velocities 
may be considerably higher. However, the near coincidence 
of red and blue-shifted gas (Fig. 8) suggests that random 
turbulent or tranverse velocities in the entrained gas domi- 
nate the measured velocity. As the line wings of both species 
are thought to be optically thin, Fig. 9 suggests a factor of 
30 variation in the relative abundance. However, there are 
other possible explanations: Girart et al. (1999) observed a 
similar factor of 10 variation in the HCO + to 12 CO J=2-l 
ratio across the NGC2071 outflow spectra. They discussed 
several possible explanations, for example changing excita- 
tion conditions, but concluded that an increase in HCO + 
abundance was the most likely cause. 

If we assume a constant excitation temperature of 30K 
and the gas is in LTE, then the results indicate Xhco+ varies 
from 2 x 10~ 8 up to 6.8 x 10~ 7 at the extreme velocities. 
This assumes the canonical value for xco — 5 x 10~° ap- 
plies throughout the gas. It represents a factor of up to 10 3 
increase in the abundance compared with the ambient gas 
in the HH2 dust clump (see previous section). If the excita- 
tion temperature is 100K, this enhancement may be as high 
as 10 4 . The optical extinction through dust mixed with the 
high-velocity gas is likely to be negligible, so external UV 
irradiation cannot easily explain the velocity-dependence of 
the HCO + enhancement. However, Taylor & Raga (1995) 
predicted HCO + abundances of up to 10 in turbulent mix- 
ing layers associated with relatively low- velocity (40kms~ ) 
shocks. This could explain both these enhancements, as well 




0.2 0.4 0.6 0.6 i 



log Vrel(km/s) 

Figure 9. Ratio of the J=l-0 HCO+ to J=3-2 12 CO line intensity 
in the region of HH2E (5 h 36 m 26.52 s , -6°47'12"), as a function 
of relative velocity, for the blue (upper histogram) and red-shifted 
(lower histogram) gas. Also shown is a fit to the average of the 
blue and red-shifted components; the slope of this line is 2.0. The 
spectra were obtained by smoothing both datasets to the same 
angular resolution (15 ") and spectral channel width (1.3b?u _1 ). 
Only the velocity regions with significant detections of both lines 
are shown. 

as the close association of the high-velocity regions with the 
low-excitation, H2-dominant knots. 

The brightness ratio in Fig. 9 shows a monotonic in- 
crease with velocity, consistent with Xhco+ x v rei- HCO + 
line wings in several outflows are found to have a relatively 
constant or even increasing brightness temperature out to 
the highest velocities (eg HH7-11 - Lizano et al., 1988; L1551 
- Rudolph, 1992; NGC2071 - Girart et al., 1999; OH231.8 - 
Sanchez-Contreras et al, 2000). By contrast, the intensity of 
CO line wings are generally found to decrease with v~ e ;' 8 (eg 
Richer et al., 2000), suggesting that the v 2 dependence of 
Xhco+ ma y be a common feature of shocked and entrained 
gas in such outflows. As depletion of HCO + occurs on a 
timescale significantly longer than the age of the flow, then 
the velocity-dependent abundance suggests either that: (1) 
a range of shock velocities exists, with the HCO + enhance- 
ment depending on the shock kinetic energy, or (2) HCO + is 
formed at the higher shock velocities, and this gas gradually 
mixes with ambient un-enhanced material. The latter might 
take place in a steady-state turbulent boundary layer (eg 
Taylor & Raga, 1995). We can use the results to estimate 
e(v), the mixing fraction in the turbulent boundary, ie the ra- 
tio of masses of initially shocked gas, m s (v), to ambient gas, 
m a (v), in this turbulent boundary. The HCO + results above 
indicate that m a (v) is approximately constant, whereas CO 
data show approximately that m s ~ v~ 2 . This would im- 
ply e(v) ~ (v/vj) 2 , where Vj is the maximum HCO + veloc- 
ity. Further modelling of such low-velocity turbulent layers 
would be of interest. 



4.3 HCO enhancement and shocks in the clump 

The relatively modest enhancement of the HCO + abun- 
dance towards the ambient velocity HH2 clump has been 
ascribed to a radiative precursor from the optically-bright 
HH knots HH2A/H (eg Wolfire & Konigl, 1993; Raga & 
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Williams, 2000). However, the current results (Fig. 7) do 
not show the limb-brightened morphology in HCO + pre- 
dicted by Raga & Williams, even though the flow axis lies 
close to the plane of the sky. Instead the results suggest a 
closer association with the low-excitation shocks, found by 
H2 and high-velocity HCO + within ~ O.Olpc of the clump. 
Luminosity arguments suggest the dust clump is heated by 
the flow itself. Furthermore, the dust morphology indicates 
that one side of the clump has been truncated by the flow im- 
pact (Fig. 2). Clearly, then, at least some parts of the clump 
have been dynamically affected by the outflow. It is perhaps 
worth noting that low-velocity HCO + clumps in the outflow 
from NGC2264G are also closely associated with low-energy 
H2 shocks and accelerated gas (Girart et al., 2000). So could 
the enhancement of HGO + at velocities close to ambient be 
due to low- velocity shocks within the dust clump, in a similar 
mechanism as evoked to explain the high- velocity enhance- 
ment? 

The HCO + linewidth in the HH2 clump is similar to 
that of the bright cores associated with VLA1 and VLA3 
(eg Choi & Zhou, 1997), even though the mass is 10-100 
times lower. We can estimate the ratio 2K/P in the HH2 
and VLA1 clumps, where K and P are the kinetic energy 
(assumed due to turbulence) and the gravitational poten- 
tial energy. Using the masses and core sizes from section 3.1 
above, and from Choi & Zhou, we find ratios of 65 and 0.6 
for HH2 and VLA1, showing that, unlike VLA1, the clump 
near HH2 is clearly not gravitationally bound. One possi- 
ble mechanism for the line broadening is disruption by the 
HH2 jet. In which case the bright ambient-velocity HCO + 
could also be caused by enhanced abundance in the low- 
est velocities of the turbulent mixing layer. Extrapolating 
the v 2 enhancement seen in the high-velocity gas down to 
v re i ~ lfcms -1 we would predict abundances a factor of 
~ 10 above normal, similar to that observed. 



5 CONCLUSIONS 

A sub-mm continuum clump is found at the end of the 
molecular outflow and 0.05pc downwind of the optically- 
bright knots in HH2. Warming of the dust appears to be 
caused by the impact of the jet, and the derived luminos- 
ity of the central peak (~ 1.2 L©) is similar to that of the 
outflow mechanical luminosity. 

The emission from HCO + can be divided into an ambi- 
ent velocity component, approximately correlated with the 
dust emission, and a high- velocity component, closely asso- 
ciated with the H2 knots. The high-velocity HCO + appears 
to form a single coherent flow linking the lower-excitation 
HH knots; this is North of the optically-bright regions, and 
we suggest it represents the main collimated jet from VLA1. 
In the highest velocity line wings, HCO + is enhanced by a 
factor of up to ~ 10 3 , compared with a factor of a few in the 
ambient- velocity gas. The enhancement is found to increase 
as v 2 el . 

Enhancement of ambient-velocity HCO + abundance is 
most prominent within 0.01 pc of the low-excitation shocks 
at the tip of the outflow jet. It is suggested that enhancement 
in this gas could be caused by the same shock mechanism 
and turbulent mixing as used to explain the high-velocity 
HCO + . Thus a UV precursor may therefore not be necessary. 
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